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Kondo Resonance in the Presence of Spin-Polarized Currents 

Yunong Qi,i Jian-Xin Zlm,2'El Shufcng Zlmng,^ and C. S. Tingi'Q 

' Texas Center for Superconductivity, University of Houston, Houston, Texas 77204 
^Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 
'^Department of Physics & Astronomy, University of Missouri, Columbia, Missouri 65211 

(Dated: February 2, 2008) 

We propose an improved method of the equation of motion approach to study the Kondo problem 
in spin-dependent non-equilibrium conditions. We find that the previously introduced additional 
renormalization for non-equiUbrium Kondo effects is not required when we use a proper decoupling 
scheme. Our improved formulation is then applied to address the spin-split Kondo peaks when a 
spin current injects into a Kondo system. 
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The Kondo effect [ij has been one of the central and 
challenging problems in condensed matter physics for 
many years. In equilibrium, the Kondo peak formed at 
the Fermi level is degenerate with respect to spin de- 
grees of freedom when a magnetic impurity is embedded 
in a non-magnetic metal. By applying a magnetic field, 
the spin degeneracy of the Kondo peak can be lifted due 
to Zeeman splitting of the impurity level. For a meso- 
scopic system such as a quantum dot, one can also use 
ferromagnetic leads to lift the spin degeneracy; this is 
because the ferromagnetic leads provide the spin depen- 
dence of the interaction (or hybridization) between the 
conduction electrons of the leads and the localized state 
in the quantum dot. To explain the above spin-dependent 
non-equilibrium phenomena, one frequently relies on the 
method of the equation motion (EOM) where an addi- 
tional renormalization in the self-energy is introduced 
without rigorous justification. 

In this Letter, we address two issues. First, we show 
that the non-equilibrium Kondo problem can be solved 
without the artificial additional renormalization as long 
as one correctly keeps the previously neglected terms in 
the EOM. We establish a new formulation for the case 
of finite and infinite U in the Anderson model. Second, 
we use the improved EOM method to predict spin-split 
Kondo peaks in the presence of a non-equilibrium spin 
accumulation. The spin accumulation has played a cru- 
cial role in the emerging field of spin electronics [2, S Ij] . 
Our prediction directly connects the separation of the 
Kondo peaks with the spin accumulation, and thus it pro- 
vides a way to determine the spin accumulation in Kondo 
systems. We emphasize an important distinction of the 
present study from previous one: the non-equilibrium 
condition generated by a voltage across the quantum dot 
and the ferromagnetic lead has a common chemical po- 
tential for spin up and down electrons, while the spin 
accumulation generates spin-split chemical potentials. 

There are a number of methods to study the Kondo 
effect in the spin-dependent non-equilibrium condition. 
The EOM approach of the Anderson model has been used 
for treating both equilibrium and non-equilibrium Kondo 



physics at low temperatures. The EOM approach in- 
cludes re-summation of low-order hopping processes and 
needs a decoupling scheme in order to obtain a closed 
analytical form. We follow the procedure introduced by 
Appelbaum, Penn, and Lacroix (APL) [5|, |6|, which is 
known to capture the right qualitative feature of physics 
at low temperatures. We note that other approaches, 
e.g., numerical renormalization group method, can also 
describe the low temperature Kondo effect [71 18| . 
The Hamiltonian of the impurity Anderson model is 



H = ^ ekacl^Cka + ^ <^dad\d„ + Ud\d„d\dg 
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Here cj^^^ and d\ are respectively the creation operators 
for conduction and d electrons at the impurity site. The 
quantities e/jo-, e^a are the conduction electron energy 
dispersion and the impurity energy, respectively. We as- 
sume that conduction electrons density of states is con- 
stant but spin dependent, i.e., pi^^a) ~ \/2D„ when 
~D^<eka<D„, U is the intra-atomic Coulomb interac- 
tion at the impurity site, and Vka represents the s-d hy- 
bridization. 

By using the standard procedure for the EOM, we ob- 
tain the impurity Green function, Gda === <C(icr|rf^3>, 
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To arrive at Eq. ^, we have made the decoupling ap- 
proximation procedure shown below. The "dynamic" av- 
erage occupation number (i.e., frequency-dependence) of 
the impurity level is defined as 

w = (n<iff -> — — f r + > -f. — 7 r , (3) 
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and the three self-energies are 
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We shah point out that the expectation value of (cqgCka), 
{dychs) which have been discarded in previous EOM 
studies are important at low temperatures since they di- 
verge logarithmically at the Fermi level as the temper- 
ature approaches to zero. Their values should be self- 
consistently evaluated through the following identities: 

{c\sCks) = / /_F_D(tj)Im<Cfcff |4<T>c^'^ , (9) 



where fpoi^) = l/[e^^'^ -I- 7] is the Fermi-Dirac distri- 
bution function, f3 ~ l/ksT, and the Green function 



<Cfc5-|c^g.> is: 



^^ , t ^^ Si,k , Fte^;ff«dff|4» 
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, (10) 



and similarly 

(4cfeff) = f{ij)Im<^Cka\4:i^dij , (11) 



with 
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Equation ([2]) is an extension to a similar result ob- 
tained by Meir et al. [9|]. Equations (2) through (10) con- 
stitute a set of the closed self-consistent equations which 
can be numerically solved. Before we carry out numeri- 
cal calculations, we point out the new ingredients in our 
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Figure 1: (Color) Spectral density pdn calculated via the EOM 
method for an infinity U Anderson model for various temper- 
atures in the absence of spin polarization. The inset displays 
the zoom-in view of the Kondo resonance near the Fermi en- 
ergy. 
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Figure 2: (Color) Spectral density pda calculated via the 
EOM method for an infinity U Anderson model for various 
temperatures for a fixed degree, P — 0.5, of spin polarization 
(a) and (b), and the Kondo resonance spfitting as a function 
of the degree of spin polarization (c) and the hybridization 
strength at a fixed P = 0.5 (d). 



formula: 1) the effective occupancy is frequency depen- 
dent, and 2) the higher order self-energy contains the in- 
termediate off-diagonal states in momentum space (e.g., 
(cjj^Cqcr)) and charge fluctuations (e.g., {d^-Cq^}). Solving 
the coupled Eq. ©-(lO) not only yields the Kondo reso- 
nance at low temperatures but also allows us to explicitly 
include the logarithmic divergence in general. 

For the case of an infinite U, the Eq. ^ takes a simpler 
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where the zero-ordered self-energy is renormahzed by a 
factor 
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Before we apply our formulation to discuss the Kondo 
resonance in the presence of the spin accumulation, we 
examine various well-studied cases by numerically solv- 
ing Eq. p^ . We choose the following parameters for 
our numerical calculation. The energy of the half-width 
of the impurity resonance in a nonmagnetic metal, Aq = 
— Im[Eocr(tt'-|-zO"'')], is taken 200 mcV unless specified oth- 
erwise. The conduction band half-width is D = IOOAq. 
The degenerate impurity level is Cd = — 6Ao. 

The simplest case is the conventional equilibrium 
Kondo problem where neither impurity state nor hy- 
bridization is spin-dependence. Figure [1] shows the impu- 
rity spectral density (pd] = Pdi) for four temperatures. 
As expected, the virtual bound level with a broad spec- 
trum and a sharp peak at the Fermi level known as the 
Kondo peak appears. The Kondo peak is suppressed 
and broadened when the temperature is increased. The 
Kondo temperature which is defined as the full-width 
of the Kondo peak is Ti<-~ exp[7re(j/2Ao]. These well- 
known results agree with many various approaches, e.g., 
the scaling analysis I 111 and the non-crossing approxima- 



tion iJl . Meir et al. [l0| and Martinek et al. [l3| has also 



used the EOM approach to derive these Kondo peaks in 
the absence of spin polarization. 

We now turn on the spin polarization. In a quan- 
tum dot, the spin polarization is introduced via the cou- 
pling to a ferromagnetic lead. In this case, one can 
parametrize the spin dependence through the hybridiza- 
tion parameter, i.e., A^ = Ao(l — trP), where the pa- 
rameter < P < 1 and (7 = ±1 for spin up (-1-1) and 
down (—1). The results arc shown in Fig. [21 When P is 
non-zero, the Kondo resonance spHts [Fig. [2l^a)-(b)], i.e., 
the peak for spin up (down) spectral density shifts below 
(above) the Fermi energy. Notice also that the intensity 
of the peak for spin up (down) is suppressed (enhanced) 
compared to the equilibrium Kondo peaks, see Fig. [2Ia)- 
(b) and Fig. [TJ These results are again consistent with 
those obtained by Martinek et al. [13[. Another way of 
introducing the spin-polarization is to apply a magnetic 
field. In this case, the Zeeman splitting of the impurity 
level becomes Co- = Cd + o'ij.bB, where B is the magnetic 
field. We find that the Kondo peaks for spin up and down 



are separated by 2/iB_B (not shown), which agrees with 
that in Ref. \M 

The agreement of our results with other approaches 
vaHdates the EOM approach of Eq. ^ or Eq. dTS]). The 
significance of our EOM approach is that it does not rely 
on the additional renormalization introduced in the pre- 
vious EOM technique [9| . The purpose of the additional 
renormalization is to account for the spin dependent level 
splitting and broadening [l0|, ll3|. The lack of rigor- 
ous justification for the existence of the additional renor- 
malization has cast a doubt for the effectiveness of the 
EOM approach for the non-equilibrium Kondo problem. 
In our improved EOM formula, we have shown that the 
correct Kondo resonance can be derived without intro- 
ducing the additional renormalization. Comparing with 
previous calculations, we have properly evaluated terms 
such as (Cq^Cka), (dgCka) througli Eq. ([5|). These terms 
make crucial contributions to the Kondo resonance split- 
ting at low temperatures. Neglecting these terms will 
lead to severe errors, which has to be recovered by the 
artificially adding an additional renormalization. 
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Figure 3: (Color) Schematic illustration of a spin-current in- 
jected from a ferromagnetic layer to a non-magnetic conductor 
containing a magnetic impurity. 



We now apply the EOM equation to calculate the 
Kondo resonance in the presence of the spin accumu- 
lation. We consider a bilayer structure where a cur- 
rent fiows from the ferromagnetic layer to the non- 
magnetic layer containing Kondo impurities as schemat- 
ically shown in Fig. [3l When a spin-polarized current 
injects into the non-magnetic conductor, a spin accu- 
mulation is built-up near the interface. Assuming that 
the ferromagnetic layer carries a spin-polarization of 
the current P, the spin accumulation is |14l . llSj 6m = 
(PjA/is/c-D) cxp(— x/A) where j is the current density, 
A is the spin diffusion length, fig is Bohr magneton, D is 
the diffusion constant, and x is the distance away from 
the interface. If we only consider the Kondo impurity 
sufficiently close to the interface, i.e., within the dis- 
tance of A, we can drop the spatial dependence of the 
spin accumulation, i.e., Sm = PjX^b l&D. To calculate 
the Kondo resonance from Eq. (13), we specify the de- 
pendence of the parameters on the spin accumulation. 
First, the spin accumulation makes the chemical poten- 
tial spin-dependent 



14l . [15[. Specifically, the chemical 



potential splitting of spin up and down conduction elec- 
5m{eDp/ p.b) = PjpX [16[, where 
Thus one should replace the Fermi 



tron is ^^ — fi 



,1 



p is the resistivity. 
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Figure 4: (Color) Spectral density pda calculated via the EOM 
method for an infinity U Anderson impurity for various tem- 
peratures in presence of spin accumulation with Apt = 2 meV. 



level in Eq. (13) by the spin-dependent chemical poten- 
tials for the spin up and down fj,^ = Ep + Pip\/2 and 
/i^- — Ep — PjpX/2, where Ep is the Fermi level. Sec- 
ond, the hybridization parameter A would also be spin 
dependent since the density of states of the conduction 
electrons is modified by the non-equilibrium electrons. 
However, the non-equilibrium electron density at very 
high current density (say 10^ A/cm ) is at least several 
orders of magnitude smaller than the equilibrium elec- 
tron density, and thus the correction to A is very small 
and we will assume that A remains spin-independent. 
Finally, the spin accumulation could lead to the spin- 
dependent energy shift of the impurity state. If one mod- 
els the interaction between the spin accumulation and 
the impurity via a phenomenological exchange coupling, 
i.e., H' = —Jex6m ■ Si, where Sm is the spin accumu- 
lation and Si is the impurity spin, the impurity level 
would be spin-split; this will be equivalent to the case 
when the impurity is subject to an effective magnetic field 
Be// = JexSrn.- Although the magnitude of the effective 
field could be respectable for a high current density, the 
local level splitting by a magnetic field has already been 
thoroughly investigated and thus we neglect the effect of 
the direct coupling between the spin accumulation and 
the impurity. Therefore, we focus on the Kondo reso- 
nances due to spin dependent Fermi levels. In Fig. 01 we 
show the spectral density p^a- for various temperatures 
with a spin-current-induced chemical potential shift. It 
is found that the spin-current influence on the Kondo 
resonance splitting is robust against temperatures. The 
more the amplitude of chemical potential relative shift 
by spin current, the more pronounced the Kondo reso- 
nance splitting between two spin channels will be. No- 
ticeably, the amplitude of the split Kondo peaks remains 
robust against the spin-induced chemical potential split- 
ting (compare Fig. 0] with Fig.[T|), which is different from 
the case of a ferromagnetic metal as discussed before (see 



Fig. [lla)-(b)). This observation is the hallmark of the 
Kondo resonance splitting from the spin accumulation 
— a purely non-equilibrium effect. 

Advanced experimental techniques such as magnetic 
tunneling into quantum dots and STM measurement 
would be able to detect the spin-split Kondo peaks. This 
will provide a high resolution detection of the spin ac- 
cumulation. For example, for a relatively small current 
density of the order of 10^ A/cm^ , the Kondo peaks will 
be separated by about SE = 0.02 meV for a typical tran- 
sition metal such as Fe (note that the width of the Kondo 
resonance, which is of the order of the Kondo tempera- 
ture). In comparison, the reliable measurement of the 
spin accumulation based on the Silsbee- Johnson tech- 
nique [l4| usually requires the current density more than 
10^ A/cm^ in order to derive sufficient large signals. 

In conclusion, we have developed an improved decou- 
pling scheme for the EOM approach in Kondo problems 
at low temperature. We show that by keeping the pre- 
viously discarded terms in the EOM approach, the ad- 
ditional renormalization is naturally implemented. We 
then apply the improved EOM approach to study the 
Kondo effect in the presence of the spin accumulation. 
We predict that the Kondo peaks are spin-split due to 
the spin dependence of the chemical potentials. 
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